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Abstract 

As a first step toward developing a reaction model that enables a comprehensive description of 
neutrino-nucleon reactions in the nucleon resonance region, we have applied for the first time a 
dynamical coupled-channels model, which successfully describes ttN, — > ttN, nN, ttttN, KA, KT, 
reactions up to W = 2 GeV, to predict the neutrino-induced meson-production reactions with 
AS = at the forward angle limit. This has been achieved by relating the divergence of the 
axial-current matrix elements at Q 2 = to the ttN — > X reaction amplitudes through the PCAC 
hypothesis. We present the contributions from each of the ttN, nN, rnrN, KA, KT, channels to the 
F% structure function at Q 2 — > limit up to W = 2 GeV. 

PACS numbers: 13.60.Le, 13.15.+g, 12.15.Ji, 13.75. Gx 
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I. INTRODUCTION 



Recent breakthrough measurements of non-zero neutrino mixing angle 9is from Daya Bay 
and RENO experiments [j], El, which are consistent with the data from T2K, MINOS and 
Double Chooz experiments [3|-|5j, indicated a possibility of the CP violation in the lepton 
sector. Now the main issue of the neutrino physics is shifting to CP phase, mass hierarchy 
as well as further precise determination of #13. For making a progress towards this direction 
by analyzing data from the next-generation long-baseline and atmospheric experiments, 
neutrino-nucleon and neutrino-nucleus scattering need to be understood within 10% or better 
accuracy, for the relevant neutrino energy region from sub GeV to a few GeV, and < Q 2 < 4 
(GeV/c) 2 [see Eq. (J3J) for the definition of Q 2 ]. This energy region covers neutrino-nucleus 
interactions of different characteristics, namely, the quasi-elastic (QE), resonant (RES), and 
deep-inelastic scatterings (DIS). Thus a combination of different expertise is necessary to 
tackle the problem. This motivates theorists and experimentalists to get together to organize 
a new collaboration, e.g., see Ref. jfij]. 

Here we are concerned with the RES region which covers the A peak and, through 
the second and third resonance regions, up to the region overlapping with the DIS region. 
Previous models for the weak single pion production off the nucleon in the RES region, some 
of them are for the A-region only, can be categorized into three kinds of approaches. Models 
of the first kind of approaches consist of a coherent sum of resonance contributions (7 10 . 
The second one additionally has non-resonant mechanisms of the tree level 

(HUTU- The 

third one considers the rescattering also so that the ttN unitarity is maintained, and such a 



model for the A-region was developed by two of the present authors [14j, [15|. These models 



for the elementary processes have been used as basic ingredients to construct neutrino- 
nucleus reaction models. Although the previous models mentioned above consider only the 
single-pion production, double-pion production is comparably important in the RES region. 
Furthermore, 77 and kaon productions also take place, and they can be a background for 
proton-decay experiments [If], [l7j]. Some models for the weak kaon productions through 
the strangeness conserving (AS = 0) reactions [l8| and the strangeness changing (AS* = 



±1) reactions [19|, |20|], belonging to the second kind of approaches discussed above, have 
been developed so far. In order to describe those meson production reactions, the reaction 
model has to take into account the coupled-channels effects and satisfy unitarity for the 
multichannel reactions. However, such a model for the neutrino-nucleon reactions has not 
been developed so far. 

In this context, our recent development of a dynamical coupled-channels (DCC) model 



is quite encouraging [211 . |22| . Our DCC model is based on a comprehensive analysis of 
ttN, syN — > ttN, tjN, KA, KT, reactions in the RES region, taking account of the coupled- 
channels unitarity including the nnN channel. An extension of the DCC model to the 
neutrino reaction is fairly straightforward. Although we need to construct a dynamical 
axial-current model for a full development, we can actually calculate the neutrino- induced 
forward (Q 2 = 0) meson production cross sections, characterized by the structure function 
F 2 , from the cross sections for ttN — > X (X = ttN, rj N, KY. . . ) by invoking the PCAC 
hypothesis. Thus, in this report, we attempt to make a first step of extending the DCC 
model to the weak sector, by calculating F 2 (Q 2 = 0) for AS = uN — > IX (I: lepton) with 
the PCAC hypothesis, thereby setting a starting point for a future full development. We 
also remark that our estimate of the magnitudes of vN — > IKY and vN — > IrjN forward 
cross sections is, for the first time, based on a model that has been rather extensively tested 
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by the data of irN and 7 A reactions in the RES region. 

The rest of this report is organized as follows: In Sec. [TTl we describe our procedure to 
calculate F 2 for the forward neutrino-induced meson production reaction using the PCAC 
hypothesis. We present numerical results in Sec. IHH followed by a summary in Sec. [TV] 



II. FORMALISM 
A. Kinematics 

First we define kinematic variables needed for the following discussions. We consider 
the inclusive l{k) + N(p) — > l'(k') + X{j>') reactions, where (I, I') = (z/ e ,e~), {p ei e + ) for the 
charged-current (CC) reactions, while (/,/') = (y e , z/ e ), (u e , u e ) for the neutral-current (NC) 
reactions. We assume that leptons are massless throughout this paper. 

In the laboratory frame, the four-momentum are defined to be 

k = (E,k), (1) 
p = (m N ,0, 0,0), (2) 
k' = (E',k'), (3) 

and p' = k + p — k'. For massless leptons, E = \k\ and E' = \k'\. The positive quantity Q 2 
is then defined by 

Q 2 = -q 2 = 4EE' sin 2 ^, (4) 

where 9 is the scattering angle of I' with respect to I, i.e., k ■ k! = cos 9; q is the momentum 
transfer between / and I', q = k! — k. Each component of the four-momentum q is denoted 
as q = (u, q) in the laboratory frame. Hereafter we call this frame FM1. 

For later use, we also define another frame, called FM2, in which X is at rest. In 
this frame, q and p are denoted as q = (u c ,q c ) and p = (En,—c[ c ), respectively, where 
E N = ^/m 2 N + \q c \ 2 and rriN is the nucleon mass. Also, we set q c = (0,0, |<f c |) so that q c 
defines the z-direction of this frame. 



B. Cross section formula of inclusive neutrino reactions at forward angle limit 



By assuming that the inclusive l{k) + N(p) — > l'(k') + X{p') reactions take place via 
one-gauge-boson exchange processes as shown in Fig. (TJ the cross sections for the inclusive 
neutrino and anti-neutrino reactions are expressed as 



dE'dQ' 



ri2 ri 
Kj pLs a 

2tt 2 



E r 



2W X a sin 2 - + W 2a cos 2 -±W 3a sin 2 t 

2 ' 2 ' m N 2 



(5) 



Here, the label a = CCu, CCu, NCu, NCu specifies the reactions; Q' is the solid angle of 
in the laboratory frame; C a = \V u d\ 2 for a = CCu, CCu and C a = 1 for a = NCz/, NCz/; 
the sign in front of Ws ya is taken to be + (— ) for v (u) induced reactions. The structure 
functions, W^ a (i = 1,2,3), are Lorentz- invariant functions of two independent variables. 
One usually chooses Q 2 and the invariant mass W = y/s = a/(p + q) 2 for the resonance 
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FIG. 1. Schematic representation of the l(k) + N(p) —> l'{k') + X(p') reactions considered in this 
paper. 



region, but chooses Bjorken scaling variable x = Q 2 / (2p ■ q) and Q 2 for the deeply inelastic 
region. In the forward limit, 9—^0, Eq. (jSJ) reduces to 



dE'dVL' 



I 2tt 2 



(6) 



The structure function Wi }Cl is expressed in terms of matrix elements of weak currents 
between the initial nucleon N and the final X, (X\J£\N). [Throughout this paper, we 
use conventions of Bjorken and Drell [23J], and any one-particle states are normalized as 
(k\k') = S(k — k').] The weak currents J a ^ are given by 



J, 



W + iV*) - (A 1 , + i A 2 ) (for a 

(Yi - iV 2 ) - (4 - iA\) (for a 

'1-2 sin 2 9 W ) V 3 - 2 sin 2 9 W V^ - A 3 (for 



a 



CCv), 
GCP), 

NO, NO). 



(7) 



Here, V* and A 1 ^ are the vector and axial currents, respectively. The superscript i — IS 
(i = 1,2,3) denotes the isoscalar current (i-th component of the isovector current). Also, 
9w is the Weinberg angle. If one evaluates (X\J£\N) in the frame FM2, then the structure 
functions are expressed as [15 



W- 



2. a 



q 2 ^ 



^d(^i^iiv)r+i(xi^iiv)i 2 )+^ 

Z He 



Q 



2 q-J 



\N) 



where we have introduced concise notation 



J2 J2( 2 ^ 3si (p+i-p"i 

V-spin X 



E 



N 



m N 



(9) 



where Ylx means summing up all possible quantum numbers and integrating over momen- 
tum f? of all final state X, and the factor 1/2 in Eq. (Q comes from taking average for the 
initial nucleon spin. 
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We now notice from Eq. (jj]) that the 9 — > limit leads to Q 2 — > 0, and thus the structure 
function for evaluating the cross section [Eq. (Q] at 9 — > reduces to 



^(Q 2 ^o) = ^^||(x|^| g .j Q |iV)P 



(10) 



Because of the vector current conservation (X | g • V a \ N) = in the isospin limit and | q c 
at Q 2 = 0, we find that 



us. 



W 2>a (Q 2 -+ 0) = i ^ • A Q |iV)| 2 . (11) 

According to Refs. [23-26|. the divergence equations for the axial currents give 

\(X(p')\q.A a \N(p)}\ 2 = flmi\(X{pi)\^\N{p))\\ (12) 

where f n (m n ) is the pion decay constant (pion mass), and if a is the normalized interpolating 
pion field. Furthermore, the matrix element (X(p')|7r a |A^(p)) at Q 2 = can be expressed as 



2Ur 



\(X(p')\fr a \N(p)}\ 2 = ^\% aN ^ x (0)\ 2 . 



ml 



(13) 



Here, T^N^xiq 2 ) is the T-matrix element of the Tc a (q) + N(p) — > X{p') reaction in the nN 
center-of-mass frame (i.e., in the frame FM2), where the incoming pion can be off- mass-shell 
q 2 ^ ml. Using Eqs. O and ([IS]), we have at Q 2 = 0, 

^ £lW)k • ^ a |iV(p))| 2 = ^(2« e ) £ |T^x(0)| 2 



2\|2 



2tt L N oo c rriN 



f 2 

J 7T 



TIUJ 



-a n a N ^. x , 



(14) 



where n^x is the total cross section of the on-shell ir a + N — > X reactions, and we have 
used the relation T^aM^xiq 2 = 0) ~ T^n^x^ 2 = m ^) ; which is a consequence from the 
PCAC hypothesis Q, and <f = u 2 . From the fact that |vr ± ) = ^(l/v^Xlvr 1 ) ± i\n 2 )) and 
1 7T°) = 1 7r 3 ) , we finally have 



r vi 



Vl 



(Tk+n^x (for ol = CCu) 



— -v-w-n^x (for a = CO) 



7Y0O 



f 2 

■ J 7T 



(15) 



—a^N^x (for a = NO, NO). 



The results so far obtained is essentially same as in the papers by Adler 28[ and also by 
Paschos et al. (29-31]. From Eqs. ()6]) and f|T5|) . one can evaluate neutrino- induced forward 
meson production reactions at 9 = using the tcN — > X total cross sections. 
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III. RESULTS AND DISCUSSIONS 



To evaluate Eq. ( TL5]) . we need inputs of ttN reaction total cross sections. In this work, 
we employ those obtained from the DCC approach developed by the authors 21, 22|. This 



approach is based on a DCC model [321 ]. within which the couplings among relevant meson- 
baryon reaction channels including the three-body irnN channel are fully taken into account, 
so that the scattering amplitudes satisfy the two-body as well as three-body unitarity. The 
scattering amplitudes of ttN — > X with X = ttN, tttcN, t]N, KA, KY< are then constructed 
through a global analysis of pion- and photon-induced ttN , rjN, KA, KT, production re- 
actions off the nucleons up to W = 2 GeV. Details of this analysis will be reported else- 
where 22| . 



We present the structure functions F 2 of the neutrino-nucleon reactions (Fig. [2] for CC 
reactions and Fig. [3] for NC reactions). Here F2 is a dimensionless quantity defined by 
F 2 = uW2- It is found that the contribution of the nN production reactions dominate 
F2 below W = 1.5 GeV, while above that energy, the contribution of tcttN production 
reactions becomes comparable with irN, indicating the importance of the irirN reactions in 
the nucleon resonance region beyond A(1232). On the other hand, contribution of r]N, KA, 
and KH reactions are much smaller [O(10 _1 )-O(10 -2 )] than that of ttN and nirN, which is 
similar to cross sections for meson production reactions with pion and photon beams. In 



the same figures, we also present results from the Sato-Lee (SL) model [14|, [15|, |33j (dotted 
curves). This model aims to describe ttN production reactions in the A(1232) region and 
thus contains only A(1232) as resonance contributions. It is noted that F2 functions for the 
SL model shown in the figures are not obtained via the PCAC hypothesis as discussed in 
Sec. Ill BL The SL model directly gives the F 2 functions because it consists of both the vector 
and axial currents, and reasonably reproduce available neutrino-induced pion production 
data [13]. Comparing dotted curves with thin solid curves that are the full ttN production 
reactions up to W = 2 GeV, we clearly see that contributions from other than A(1232), e.g., 
higher resonances and/or backgrounds, become relevant above W = 1.3 GeV. Also, the good 
agreement between the thin solid and dotted curves for W < 1.3 GeV indicates a reliability 
of calculating F2 from the ttN — > X total cross sections with the PCAC hypothesis. 

It is noted that, even above the A(1232) region, the F 2 function still has bump structures 
and is not a monotonous function in W . This non-monotonic behavior of F 2 comes from high- 
mass nucleon resonances, which are expected to exist up to W ~ 2.5 GeV. An appropriate 
treatment of such behavior, as successfully done in our DCC approach, may be crucial for 
reducing systematic errors in atmospheric and accelerator experiments to determine neutrino 
parameters. 



IV. SUMMARY AND PROSPECTS FOR FUTURE DIRECTIONS 

As a first step toward developing a reaction model that enables a comprehensive de- 
scription of neutrino-nucleon reactions in the nucleon resonance region, we have ap- 



regi 

plied for the first time the DCC approach developed in Refs. [21|, |32| to the neutrino- 
induced forward meson-production reactions off the nucleons, I + N — > I' + X with 
(/,/') = [y e , e~), {p e , e + ), [y e , z/ e ), [p e , z/ e ) and X = ttN, htcN, r/N, KA, KT,, in the energy 
region from the irN threshold up to W = 2 GeV. This has been achieved by relating diver- 
gence of the axial-current matrix elements (X|c^j^|iV) at Q 2 — > to the ttN — > X reaction 
amplitudes from the dynamical coupled-channels model through the PCAC hypothesis. 
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CC v-proton / CC v-neutron CC v-neutron / CC v-proton 




W (GeV) W (GeV) 



FIG. 2. (color online) VF-dependence of the Fi function for CC neutrino-nucleon meson-production 
reactions at the limit Q 2 — > 0, plotted for W from the l7r-production threshold up to 2 GeV. The 
left (right) panel is for the CC ^-proton or CC z/-neutron (CC ^-neutron or CC i>-proton) reactions. 
Each curve is: (thick solid curves) Total contribution from X = ttN, ttttN, rjN, KA, KT, production 
reactions; (thin solid curves) contribution from X = ttN only; (dashed curves) contribution from 
X = ttttN only; (dashed-dotted curves) contribution from X = r/N only; (dashed-two-dotted 
curves) contribution from X = KA only; (two-dashed-dotted curves) contribution from X = KT, 
only. As a comparison, results from the SL model in which F 2 is directly calculated without 
relying on the PCAC hypothesis, are also shown as dotted curves. 



We have presented the F 2 structure functions for I + N — > I' + X and investigated 
contributions of production reactions for each X. It is found that above W = 1.5 GeV, the 
contribution of 7T7rA" production reactions becomes comparable with ttN. Also, our results 
suggest that a naive extrapolation of the DIS cross sections down to the nucleon resonance 
region, which is often performed in analyses of atmospheric and accelerator experiments, may 
be better to be replaced by more realistic reaction cross sections for precise determination 
of neutrino parameters. 

The next step will be extending our dynamical coupled-channels model to directly analyze 
neutrino reactions without relying on the PCAC hypothesis, so that we can investigate 
neutrino reactions at any finite Q 2 . This project is underway and will be reported elsewhere. 
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NC v-nucleon / NC v-nucleon 
10 i — i ■ 1 ■ 1 ■ 1 ■ — i 




FIG. 3. (color online) PE-dependence of the F2 function for NC neutrino- nucleon meson-production 
reactions at the limit Q 2 -4 0, plotted for W from the l7r-production threshold up to 2 GeV. The 
meaning of each curve is same as Fig. [2j 
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